Interaction-mediated surface state instability in disordered three-dimensional 
topological superconductors with spin SU(2) symmetry 
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We show that arbitrarily weak interparticle interactions destabilize the surface states of 3D topo- 
logical superconductors with spin 5(7(2) invariance (symmetry class CI), in the presence of non- 
magnetic disorder. The conduit for the instability is disorder-induced wavefunction multifractality. 
We argue that time-reversal symmetry breaks spontaneously at the surface, so that topologically- 
protected states do not exist for this class. The interaction-stabilized surface phase is expected to 
exhibit ferromagnetic order, or to reside in an insulating plateau of the spin quantum Hall effect. 
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The existence of novel delocalized surface states is a 
key signature of 3D topological phases of matter [lJ-Q - 
These states envelop a fully-gapped, yet topologically 
"twisted" bulk and can display exceptional properties 
such as the quantized magnetoelectric effect and Majo- 
rana fcrmions A complete classification 0, Q for 
(effectively) non-interacting particles has demonstrated 
that only five classes of topological phases and associ- 
ated surface states arise in 3D. 

An important development Q has been the incorpo- 
ration of disorder effects on 2D surface states. This 
is crucial because the terminating facets of a bulk 3D 
crystal inevitably host structural defects and impurities. 
The topologically nontrivial bulk is linked (H to an effec- 
tive low-energy surface theory of 2D Dirac fermions, per- 
turbed by random impurity potentials 0, H| • Each of the 
five classes of 3D topological phases is "protected" from 
the effects of time- reversal invariant (i.e., non-magnetic) 
disorder, in the sense that at least one surface Dirac wave- 
function escapes Anderson localization 0, [B| . 

Unlike uniform plane waves, the extended 2D energy 
eigcnstates enveloping a surface-disordered topological 
phase exhibit wild spatial amplitude fluctuations. These 
arise from quantum interference due to multiple impu- 
rity scattering, and manifest in the local density of states 
(LDOS). The pattern of LDOS fluctuations presents an 
intricate structure, characterized by an infinite set of scal- 
ing dimensions associated to interwoven fractal measures 
of the surface, a feature known as multifractality The 
evasion of localization in favor of multifractal scaling is 
rare in 2D, and is a signature of topological protection in 
the presence of disorder [g[ . 

In this Letter, we show that topological protection 
can be undermined by interparticle interactions. In par- 
ticular, we study the combined effects of multifractal 
LDOS fluctuations and interactions upon the surface An- 
dreev bound states of 3D topological superconductors. 
Because the bulk condensate screens the long-ranged 
Coulomb force, surface quasiparticles interact only via 



short-ranged potentials. In the clean limit, the vanish- 
ing density of states for the 2D Dirac surface band im- 
plies that weak short-ranged interactions are irrelevant, 
i.e. the surface states remain "protected." However, it is 
known that disorder-induced LDOS multifractality can 
amplify interaction effects, such as pairing correlations 
near the superconductor- insulator transition Q. With 
physics dominated by its surface, the complete picture of 
a 3D topological phase must incorporate both disorder- 
induced quantum interference and interactions 0]. 

Specifically, we demonstrate that arbitrarily weak in- 
teractions (consistent with bulk symmetries) destabi- 
lize the non-interacting surface states of 3D topologi- 
cal superconductors with spin SU(2) symmetry (class CI 
El), in the presence of non- magnetic disorder. Multi- 
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FIG. 1: Phase portrait sketch for the surface physics of a 
3D time-reversal invariant, spin 5(7(2) symmetric topological 
superconductor. The vertical axis is the interaction strength 
?7 [Eq. pop ], while the horizontal axis measures non-magnetic 
disorder. Although the non-interacting system has a disorder- 
stabilized phase with delocalized ("protected") surface states 
(II), it is destroyed by arbitrarily weak interactions [Eq. (fTTj) ]. 
Instead, at zero temperature, we expect that the surface ex- 
hibits broken spin symmetry [U > (III)], or the spin 
quantum Hall effect [U < (IV)]. In either scenario, in- 
teractions break time-reversal symmetry spontaneously. 
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fractal LDOS fluctuations enhance the interactions, fa- 
cilitating the instability. We argue that time-reversal 
symmetry breaks spontaneously at the surface, so that 
"protected" surface states do not exist in this class. De- 
pending upon the sign of the relevant interaction coupling 
U [see Fig. [TJ Eqs. (HUJ) and ([H])]. the surface should de- 
velop ferromagnetic order, or enter an insulating plateau 
state of the spin quantum Hall effect [i| . Our result pro- 
vides impetus to identify a suitable material for the class 
CI bulk as an avenue to realize the spin quantum Hall 
phase. A similar analysis for the 3D topological super- 
conductor class AIII will be published elsewhere [Ioj |. 

Three of the five 3D topological symmetry classes can 
be realized as time-reversal invariant superconductors, 
distinguished by the degree of electronic spin conserva- 
tion. In a 3D topological superconductor, Cooper pairing 
leads to a fully-gapped quasiparticle band in the bulk, 
associated to an integer- valued winding number v 0, Q . 
The modulus \v\ equals the number of flavors (or "val- 
leys" ) of 2D quasiparticle bands that appear at the sam- 
ple surface, with energies that infiltrate the bulk gap. 
In the clean limit, the surface states exhibit a massless 
Dirac character at low energies; the Dirac point appears 
precisely at the bulk chemical potential (inside the gap) 
due to particle-hole symmetry. 

In this paper, we study a universal low-energy model 
for the 2D surface states of a 3D class CI topological 
superconductor. In contrast to the spin-orbit-coupled 
Z2 topological insulators, a CI superconductor has full 
spin SU (2) symmetry. The non-trivial topology arises 
through the entwining of orbital degrees of freedom, in- 
cluding non- (simple) s-wave pairing (III . 12 1. For class 



spin- 1/2 electrons in the bulk translates into the follow- 
ing chiral condition on the surface 
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CI v is even because Dirac surface bands appear in time- 
reversal conjugate pairs 0, [l2|, EH- We consider the 
generic case with \v\ — 2k, k £ {1,2,3,...}. Neglect- 
ing interactions, the Hamiltonian for the surface theory 
is 



H D = / d 2 r 0t {-a- ■ [iV - A,(r)ty } 0. 



(1) 



The fermion field ip is a complex Dirac spinor with pseu- 
dospin er G {1,2} and valley v G {l,...,2fc} indices, 
i.e. tjj — > ip a ,v when all indices are displayed. The pseu- 
dospin components ipi. v and ^2. v are linear combinations 
of the Nambu elements c^ tV ^\ and ct - x ,. These annihi- 
late (create) spin up (down) electrons in valley v (v). 
(Under time-reversal, v and v interchange.) The indices 
{A, A'} label additional orbital (e.g. sublattice) degrees of 
freedom, whose precise interpretation depends upon bulk 
microscopies. A 3D class CI lattice model with v = ±2 
appeared in Ref. (l2l ]. 

For a 3D topological superconductor, a key conse- 
quence of the non-trivial bulk is the special form that 
time-reversal symmetry adopts on the surface. If we 
write Ho = ip^hip, with h the single-particle Hamilto- 
nian operator, then the usual time reversal symmetry for 



(2) 



This implies that any surface disorder that does not 
break time-reversal (including non-magnetic impurities) 
can manifest only as a random vector potential in the 
low-energy Dirac description. [Recall that ip in Eq. (TTJ) 
carries U(l) spin, rather than electric charge. In this lan- 
guage, vector potentials couple to time-reversal invariant 
spin and valley currents.] A homogeneous perturbation 
such as a chemical potential shift, or a time-reversal in- 
variant pairing of the surface quasiparticles can be elim- 
inated by a gauge transformation. Moreover, an energy 
gap (Dirac mass term) cannot appear at the surface of 
a topological superconductor unless time-reversal is bro- 
ken. For class CI, non-magnetic disorder induces scatter- 
ing between the 2k valleys, in the form of the non-abelian 
valley vector potential Ai(r)t^ in Eq. (JT|). Here i\ de- 
notes a,2kx2k generator of the group Sp(2k) . (The group 
is symplectic due to the spin symmetry |l3(.) In the ab- 
sence of interactions, elastic scattering due to vector po- 
tential disorder begets delocalized, multifractal surface 
states, many properties of which can be computed ex- 
actly via conformal field theory (CFT) H0[l5j]. 

We first consider the effects of disorder upon the non- 
interacting surface states. Below we describe the physics 
and main idea of the CFT method. A technical summary 
can be found in Ref. [l3^ . while a more comprehensive 
discussion will appear elsewhere (lp| . The spatial char- 
acter of the surface state wavefunctions (localized ver- 
sus extended) can be ascertained via disorder-averaged 
moments of physical observables, such as the conduc- 
tance or the local density of states (LDOS). To facilitate 
this, we replicate ip a .v — ^ Vv.w.a, where the replica index 
a G {1, . . . , n}, and we are to take n — > at the end of the 
calculation [5|. Symmetry is the primary tool employed 
in the following, so we will rewrite Eq. ([1]) in a manifestly 
symmetric form. We decompose ip and ift into "left" C 
and "right" K fields, 

{TZ^,v,a,Ti-l,v,a} = \^2,v.a, tp\^> <a (k 2 )v> ,v } ■ (3) 

Here and below, repeated indices are summed. £ StV . a de- 
notes a 4n£;-component spinor; the index s (v) transforms 
in the fundamental representation of the spin SU(2) [val- 
ley Sp(2k)] symmetry. We also define 



L = C 



R = K T is 2 k 2 



(4) 



\-a' v and Rq' v transform in the conjugate representations 
of the spin and valley symmetry groups; s 2 and k 2 
are spin and valley antisymmetric Pauli matrices [lij ]. 
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Eq. (UJ) can be rewritten as Hd = Hq + SHd, where 

H =i [ d 2 r [LdC+RdK] , (5) 

(6) 



SH D = Jd 2 r [j^ + J^] 



In Eq. (JS|), we have switched to complex spatial coordi- 
nates {z, z} = x±iy, {d, 8} = \{d x ^id y ). The valley dis- 
order appears in Eq. (|6]), where {Ai, Ai} = —i(Af^fiA^). 
This couples to the valley Sp{2k) current, which has the 
holomorphic component J* = — (i/2)LV K C. 

Eq. ([5]) is manifestly invariant under chiral (indepen- 
dent left and right) spin SU (2), valley Sp(2k), and replica 
SO(n) transformations. The symmetry group enlarges to 
SO(4nk) if we include operations that mix all three in- 
dex types. This free fermion theory is equivalent to the 
SO{Ank) 1 Kac-Moody CFT 0. The latter has a spe- 
cial property known as a conformal embedding rule 17- 
Hij , which gives a decomposition into a "product" of two 
other CFTs: Sp(2n)k, associated to the (spin) x (replica) 
invariance of Eq. (|5]l. and Sp(2k) n , associated to the val- 
ley symmetry. 

The derealization physics of the non-interacting sur- 
face with Hamiltonian Hq + SHd is the same as in 
Refs. 14, 15| . which dealt with 2D Dirac fermions cou- 
pled to a random SU(N) vector potential. Disorder is 
a relevant perturbation to the clean fermion theory flij ]. 
Crucially, the impurity potential couples only to the val- 
ley current J* in Eq. ©. This leads to a "fractionaliza- 
tion" of the original SO{4nk) 1 CFT: the valley Sp{2k) n 
sector localizes, leaving behind the "critical" (delocal- 
ized) spin-replica Sp(2n) k sector [l4j|- The latter is used 
to compute the scaling behavior disorder-averaged ob- 
scrvablcs. Even in the absence of interactions, disorder 
localizes all surface states away from zero energy [20| ; 
this is different from the case of a single Dirac fermion 
on the surface of a 3D topological insulator 21 1. How- 
ever, the localization length diverges upon approaching 
the chemical potential, so that the zero energy state at 
the Dirac point remains completely delocalized ("topo- 
logically protected" ) . 

The disorder-induced spatial fluctuations of the LDOS 
^(e.r) are encoded in the multifractal spectrum r{q) 
[a, Iff. The r(q) spectrum measures the sensitivity of ex- 
tended wavefunctions to the sample boundary. A large 
L x L area of the surface is finely partitioned into a grid 
of boxes of size a <C L. One then defines the box proba- 
bility fi n and inverse participation ratio V q , 



fx n (e) = J d 2 rv(e,r), V q (e) = ^ 



fJ-n(s) 



(7) 



An 



where A n denotes the n th box and v = ^ [ii is the global 
DOS. When e is tuned to a critical derealization energy 
(such as a mobility edge), V q ~ (a/L) T< - q \ where the ex- 
ponent r(g) is both self-averaging and universal [22[. The 



multifractal spectrum thus provides a unique fingerprint 
for spatial fluctuations in a particular symmetry class. 
In the field-theoretic description, the q th moment of the 
disorder-averaged LDOS (q G {1,2,3...}) is associated 
to a particular composite operator O q , with scaling di- 
mension A q . The set of such dimensions determines the 
multifractal spectrum via r(q) = 2(q— 1)+A q — qAi 
By contrast, localized states are insensitive to the sample 
boundary for sufficiently large L and have t(<j) = 0. 

For the class CI surface, we have identified the oper- 
ators that represent disorder-averaged LDOS moments; 
these are a subset of the primary fields in the Sp(2n)k 
CFT. As a result, we obtain the exact disorder-averaged 
multifractal spectrum at zero energy [ToL 13 1, 



r{q) = («-!) 



2(k + 1) 



(8) 



For k = 1, Eq. © agrees with previous calculations 15|; 
the form for general k is new. One of the main results of 
this paper, Eq. ([8]) proves that the non-interacting surface 
states at the bulk chemical potential remain delocalized, 
a consequence of the bulk topological order. 

Now we turn to interparticle interactions. Robust sur- 
face states must be protected from the combined ef- 
fects of both disorder and interactions. In a weakly- 
interacting fermion gas, the low-energy behavior of the 
density of states completely determines the importance of 
short-ranged interactions. The lowest-order (tree level) 
renormalization group (RG) equation for a generic four- 
fermion coupling U is [10( 



dlnU/dl = Ai — A 



(U) 



0(U). 



(9) 



where I denotes the log of the RG length scale such as the 
system size. In a clean 2D system, A^ = 2Ai, with Ai 
the scaling dimension of the LDOS. For the clean Dirac 
surface band, Ai = 1, so that weak short-ranged interac- 
tions are strongly irrelevant. By contrast, a negative Ai 
(due, e.g., to a van Hove singularity) would imply that 
U is relevant, signaling a potential instability. With im- 
purities present, the exponents Ai and A^ denote scal- 
ing dimensions of the disorder-averaged LDOS and four- 
fermion interaction, respectively. The latter satisfies the 
lower bound A^ > A2 [Io| . where A2 is the dimension 
of the second LDOS moment that determines r(2). The 
crucial point is that A2 is independent of, and strictly less 
than 2Ai for a multifractal delocalized state in a disor- 
dered system [23j ■ Impurity-mediated LDOS fluctuations 
can therefore amplify short-ranged interaction effects, by 
increasing the overlap of single particle wavefunctions in 
local regions. This is particularly relevant for an interac- 
tion U that saturates the bound A^* 7 ' = A 2 < 2Ai. 

Physically, we expect that the important interactions 
include a spin exchange channel (because spin is a con- 
served hydrodynamic mode) and a Cooper pairing in- 
teraction (because disorder respects time-reversal). The 
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former is written —S-S, where S denotes the spin den- 
sity. As discussed below Eq. ((2]), pairing of the surface 
quasiparticles does not open a gap unless time-reversal 
is simultaneously broken. The latter occurs when the 
Dirac mass operator m = ip^a 3 ^ develops an expecta- 
tion value. This can be understood explicitly in the 3D 
CI topological superconductor lattice model of Ref. [l2| , 
which features real d-wave pairing in the bulk. In that 
model, m is interpreted as a sum of pairing operators: 
in ~ — icjcj + ic^Cf, where c s annihilates a lattice elec- 
tron. Crucially, in is odd under time-reversal, due to 
the factors of i. Thus, a non-zero expectation (m) ^ 
would imply ("d + is") pairing of the surface quasipar- 
ticles, opening an energy gap and breaking time-reversal 
symmetry. The resulting state is an insulating plateau of 
the spin quantum Hall effect (see below). An attractive 
Cooper pairing interaction can be written as — m 2 . 

To keep the analysis general, we enumerate all four- 
fermion interactions that preserve the bulk symmetries 
[time- reversal invariance, spin SU(2), and valley Sp(2k) 
symmetry] . This necessitates the incorporation of a third 
interaction channel JgJg-, where is the holomorphic 
spin current. The replicated interaction Hamiltonian for 
the CI surface is UM 



U{m a m a -AS a -S a )+VJ 1 Sa J 1 Sa 
1 



+W ( 3m a m a + 4S a ■ S a - ^J2 a 



(10) 



The interaction strengths {U, V, W} arc defined so as to 
couple to RG eigenoperators, in the presence of disorder. 
In the minimal two valley realization (k = 1), the W- 
channel interaction does not exist. For that case only, 

Jsa^Sa = 3m a m a + 4S a ■ S a . 

Our task is to evaluate Eq. (j9j in the disordered, non- 
interacting CI surface theory for the three interaction 
operators in Eq. (JTDJl 0. Using the Sp{2n) k CFT, we 
have found that one particular operator controls the scal- 
ing of both the second LDOS moment and the interaction 



U, leading to A y 2 



(u) 



A 2 = 0, while Ai = l/2(& + 1) 



Tol 13j . The main result of this paper follows, 



dU_ 



u 



2(k + l) 



o (u 2 ) , 



(11) 



which implies that the interaction U in Eq. (|T0)) grows 
at longer wavelengths, destabilizing the non-interacting, 
dirty surface, for any number of 2k valleys. By contrast, 
the other interactions V and W remain irrelevant for any 
k, satisfying = ^ = UM- 



dl 2(fe+l) ' dl 20+1) 

We conclude that while weak interactions are suppressed 
in the clean limit by the vanishing density of states at 
the Dirac point, surface disorder strongly renormalizes 
the interaction channel U, making it relevant. 

Eq. (|11[) can be understood as an enhancement of inter- 
action matrix elements in the eigenbasis of the disordered 



theory: local accumulations of the DOS due to wavefunc- 
tion multifractality induce stronger interactions between 
the surface quasiparticles. The amplification of the par- 
ticular interaction channel U over the others signals the 
instability of the non-interacting surface to spontaneous 
time-reversal symmetry breaking. From Eq. (|10[) , we an- 
ticipate (at least local) ferromagnetic order (S) ^ when 
U — > +oo. Without time-reversal symmetry, the surface 
is not "topologically protected" [lrl3| , and we expect An- 
derson localization of all surface states (HQ. However, 
we cannot rule out an exotic metallic phase when spin 
symmetry is also broken [25|. By contrast, U — > — oo 
should cause Cooper pairing of the surface quasiparticles. 
Treating the relevant interaction in mean field theory, 
one replaces m 2 — >• 2{m)ip^ a 3 ip in Eq. ([TO)) . A non-zero 
Dirac mass opens an energy gap, producing an insulat- 
ing surface. Time-reversal symmetry is broken because 
(in) implies surface pairing at a non-zero supcrfluid 
phase angle with respect to the bulk. 

To lowest order in (1/fc), Eq. (TTTJ agrees with a per- 
turbative result (26j obtained using the non-linear sigma 
model 0, [l3|. The calculations in Ref. [26[ were per- 
formed in the context of a non-topological 2D system 
of gapless superconductor quasiparticles, subject to dis- 
order and interactions with spin SU (2) symmetry and 
time-reversal invariance. The Sp(2n)k CFT employed 
here has a sigma model representation with the same 
structure, but augmented with a Wess-Zumino-Witten 
(WZW) term In the k 1 limit, this model be- 

comes weakly-coupled, and the WZW term can be ig- 
nored. The results of Ref. [2(| therefore provide a non- 
trivial check of our analysis in the many- valley limit. In 
addition, at one loop in the sigma model calculation, 
RG flow equations beyond linear order in the interac- 
tion strengths can be obtained, because the sigma model 
treats interactions non-perturbatively via RPA and BCS- 
type summations. For the 2D class CI quasiparticle sys- 
tem, the sigma model generically predicts an instabil- 
ity of the "metallic" phase signaled by the divergence of 
the spin exchange or BCS pairing interaction strengths 
This provides evidence for the absence of an 
interacting, time- reversal invariant fixed point. 

The insulating state that occurs for (m) ^ preserves 
spin SU(2) symmetry. This state resides in a plateau 
of the so-called spin quantum Hall effect @, analogous 
to the "half-integer" quantum Hall phase at the sur- 
face of a 3D Z2 topological insulator with broken time- 
reversal symmetry [J Il2| . The quantized spin Hall con- 
ductance [H is a% y = \ (| ) 2 p, with p = fcsgn(m) if val- 
ley symmetry is unbroken on average (i.e., after disorder- 
averaging). If valley symmetry remains broken even after 
disorder-averaging, then p G { — k, —k + 2, . . . , k — 2, k}; 
see also Ref. [9( . Our results are summarized in Fig. [TJ 

In conclusion, we have demonstrated that interactions 
destabilize class CI disordered surface states in 3D. We 



5 



have argued that time-reversal symmetry breaks spon- 
taneously, and that the CI topological superconductor 
surface enters into either a ferromagnetic or a spin quan- 
tum Hall phase. These are expected to be interaction- 
stabilized Anderson insulators. The other 3D topological 
superconductor classes AIII and Dili also admit WZW 
CFT descriptions Q. The minimal surface state (single 
Dirac valley) realization for each of these is stable against 
disorder and short-ranged interactions Results for 

class AIII with multiple valleys will appear elsewhere , 
while class Dili is an important topic for future work. 
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SUPPLEMENTARY INFORMATION 

In this supplement, we provide the minimal details necessary to reproduce the results in the main text; a com- 
prehensive discussion and analysis of class CI and AIII topological superconductors will appear elsewhere We 
first summarize the symmetry structure of the class CI surface state theory. Then we provide the derivation of the 
multifractal spectrum and interaction scaling dimensions. We close with a review of the many-valley limit and the 
relation to the non-linear sigma model. 

SYMMETRIES OF THE SURFACE STATE THEORY 

In Eq. (1) of the main text, the Dirac fermion ^^(r) is a surface-confined projection °f a bulk field v I' crjt ,. T (r, z) 
that lowers the spin angular momentum by one. Relative to the surface, the bulk fermion carries an additional index 
r e {1, 2}. The four a ® r components of ty v are linear combinations of electron operators c^ yVt \ and ct - x , , where 
(v,v) denote a pair of valleys related by time-reversal, and {A, A'} run over additional orbital (e.g. sublattice) labels 



In the bulk, a physical time-reversal operation sends i — > — i (complex conjugation) and 

&t,v ->■ —Mc± t v, c^ v ->■ Mc tjS . 

In this equation, M — > M\^\i is a symmetric, unitary matrix in orbital labels. Time-reversal induces the transformation 
^ — > — i<r 3 (2P — 1)^', where P = (1 + n • r)/2 projects onto a certain r-spin direction n. (The latter is basis-dependent 
and determined by microscopies At the surface, the r-spin becomes "locked," with ip(r) ~ P^r, z = 0) [!, H|. 
This type of projection always occurs at the surface of a d-dimcnsional topological insulator or superconductor, 
producing an anomalous state that is "half" of a normal (d — l)-dimensional system. 

For the surface theory in Eq. (1), time-reversal therefore appears as the antiunitary transformation 

V; -> -ia 3 ty t ] T , i^-i. (SI) 

Spin SU(2) symmetry requires invariance under U(l) if) — > e l9 ip and particle-hole ip — > ia 1 ^ 2 ["0^] T transformations; 
the latter encodes a 7r rotation by S x . Imposing these upon Eq. (1) restricts the disorder to the Sp(2k) valley 

vector potential A 2 tJ., where — k 2 [t^,] k 2 = t l K . In these equations, k 2 denotes the valley Sp(2k) invariant tensor 
(antisymmetric 2k x 2k block Pauli matrix). 

CONFORMAL EMBEDDING 

The disordered, non-interacting class CI surface state theory is described by the Sp(2n)k Kac-Moody CFT. To see 
this, we note that the free fermion Hamiltonian in Eq. (5) of the main text is equivalent to the SO(4nk)i current 
algebra @. (Since we are discussing non-interacting fermions at this stage, we can trade the 2D Hamiltonian for 
a 2+0-D Grassmann field action.) The conformal embedding SO(4nk)i D Sp(2n)k © Sp(2k) n implies that the 
associated stress tensor T(z) has the Sugawara decomposition 0] 

T (z) = 3r^+iy [: « :(*)+: J% R J% R = («)] , (S2) 

where 

J* = I ^] ^Kal = _ L a t^.£ a , Jg R = — ^Ltg R C. (S3) 



The rcplica-summcd valley current J % K (z) satisfies the Sp(2k) n algebra. In Eqs. (|S2[) and (|S3[) . Jg R generates Sp(2n) 
spin x replica space transformations; (t^V 5 aa> is a suitable 2n x 2n matrix. The current Jsr{ z ) satisfies the Sp(2n)k 
algebra. 

The impurity potential {Ai, Ai\ in Eq. (6) of the text couples only to the valley Kac-Moody current {J l K , J*}. As 
the disorder renormalizes towards strong coupling, it localizes the valley Sp(2k) n sector. To obtain Eqs. (8) and 
(11), we express the associated operators in terms of Sp(2n)k primary fields and descendants. The most relevant 
components govern the leading scaling behavior. 
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MULTIFRACTAL SPECTRUM 

We consider first the r(q) spectrum at level k — 1. Primary fields are labeled by the Sp{2n) fundamental weights 
{A g }, q E {l,...,n} @. Weight A g corresponds to a rank g, fully antisymmetric tensor fi^'^ ...^ p satisfying the 
traceless condition (s 2 ) 1 - 7 ^^...^] = 0. Here, i = {s,a} is a product of spin s and replica a indices; (s 2 )^ is the Sp(2n) 
invariant tensor. In the holomorphic sector, the q th LDOS moment must correspond to a tensor with q distinct replica 
indices. We therefore associate 



^ n [{1 1 a}{ S2 ,2}...{ S ,,g}]( Z )^[{]' 1 ,l}{4^}---{<,?}]( 2; ) 

x ( Z s 2 s 3 ) Sl ' s 'i(is 2 s 3 ) S2 ' s 2 x • • • x (is 2 s 3 ) s "<. (S4) 

The left-hand side is the disorder- averaged g th moment. The right-hand side is a diagonal primary field labeled by 
Ag. Here we have written each Sp(2n) index as the spin-replica product {s, a}; ip'tp «S 3 is the spin-projected LDOS 
[8f for the surface state. We find that the most relevant contribution to the q th LDOS moment is associated to the 
same Sp(2n) representation A 9 , for any k > 1. The scaling dimension of the q th moment is then given by 0] 

A q = (2q-q 2 )/2(k + l) (S5) 

in the replica limit | n->0. Eq. (8) follows. 

INTERACTION DIMENSIONS 

The flow equations for each of the three interaction couplings U, V, and W in Eq. (10) of the main text appear as 
in Eq. (9), with Ai given as above and A2 ^ the scaling dimension of the associated field. We rewrite Eq. (10) 
exploiting Sp(2n) [10( and SU(2) Fierz identities: 



Hj=J d 2 r [2u.r Ka .r Ka + v.r s Ji a + 2wr v ?J^] , (se) 



where Jg a = — L a s 7 £ a denotes a replica-resolved spin current (no sum on a is implied). Here s 7 denotes a spin space 
Pauli matrix. The holomorphic half of the last term in Eq. (|S6j) is I^? a = L^P £ v > a + K'^^Sa- ^ we assume that 
the disorder-averaged theory is invariant under both spin SU(2) and valley Sp(2k) transformations, then the three 
channels U, V , and W in Eq. (|S6|) exhaust the possibilities for four-fermion interactions. 



As V couples to a KM current-current perturbation, A^ = 2 0. The U interaction involves the valley current 
J* a . This is not a valley KM current, but rather a product of Sp(2n)k and Sp(2k) n primary fields. [Replica-resolved 
components cannot be extracted from the KM current J l K in Eq. (|S3|) .] In the Sp(2n)k theory, J* a corresponds to 
a second rank tensor with equal replica indices, antisymmetrized over spin (to obtain a singlet). The only choice is 

(2) 

J l Ka {z) <=> ^[ { ; a}{s , a}] (z), i.e. the same representation A2 that determines the scaling of the second LDOS moment 
[Eqs. (|S4| and (fS5|t with q = 2]. As a result, A2 = A2 = 0. Finally, the W interaction operator I^7 a is a second rank 

tensor field with equal replica indices, symmetric in spin. The only choice has weight 2Ai, leading to A^' 1 = 2/(fc + 1) 
Via Eq. (9), we obtain Eq. (11) and the equations for V and W quoted in the text. 

MANY VALLEY (LARGE k) LIMIT: REVIEW OF THE CLASS CI FINKEL'STEIN NLctM RESULTS 



In Ref. [ill ], the effects of interactions upon gapless quasiparticles in disordered, non-topological 2D superconductors 
were considered using the Finkel'stein non- linear sigma model framework [I2I ]. In particular, the author studied 
symmetry class CI, appropriate to a spin singlet superconductor possessing spin SU{2) symmetry and time-reversal 
invariance in every realization of the disorder. The 2D topological surface state Sp{2n)k CFT studied in the present 
paper also possesses a sigma model description with the same structure, but augmented by a WZW term 0, SSLl3l ■ In 
the limit of many valleys i > 1, the WZW term can be ignored in the first approximation, yielding a "metallic phase" 
with a large spin conductance proportional to k, independent of the disorder (9)]. The advantage of working directly 



8 



in the sigma model framework is that interactions are treated to all orders, via RPA and BCS-type rcsummations. 
The disadvantage is that the sigma model becomes strongly coupled for small k G {1,2,3, . . .}, so that the WZW 
term cannot be ignored and perturbation theory becomes unreliable. 

The only interactions that appear in the CI sigma model treated in Ref. are short-ranged, and include a spin 
triplet exchange coupling j t and a Cooper pairing interaction j c . Here, we define jt,c > for repulsive interactions in 
each channel; these are dimensionless couplings in the limit of vanishing disorder strength. Other (e.g. electric charge 
density-density) interaction channels do not couple to conserved hydrodynamic diffusion modes, and are strongly 
irrelevant in the sigma model framework. 



The one-loop RG equations for 7 tjC are [11 



^ = -^7c(l-7t)(l-27t), 

d A (S7) 

= 2 {~ 3 ^ - 2 ^ + 3 7c [log(l - It) + It]} ~ ll 

In these equations, A is proportional to the dimensionless inverse spin conductance. For weak interactions, these 
equations can be linearized in the coupling strengths. The result is 

rfln> _ A d\ny> _ 3 A 
dl ~ 2' dl ~ 2' 

where 7 r = 7 C — 3jt, and j, = 7 C + 7t- Since A oc 1/k 0, 0], we find that Eq. JS8| matches the RG equations for the 
relevant {U O 7 r ) and irrelevant (W f-> 7) coupling strengths in the WZW surface state theory studied in the present 
paper, Eq. (11) and the text following, valid to lowest order in 1/k. The relevant sigma model coupling j r is indeed 
the difference of repulsive singlet pairing and triplet exchange interactions, just as we found in the WZW model. 

Integrating the full flow Eqs. (|ST[) numerically, one sees that the metallic phase is generically destroyed by one 
of two instabilities: either the triplet exchange interaction flows to minus infinity 74 — > —00 (suggesting Stoner 
ferromagnetism) , or the Cooper pairing strength flows to minus infinity 7 C — > — 00, indicating a residual pairing 
instability for the gapless quasiparticles. We do not discuss here the back-reaction of the interactions upon the spin 
conductance [ll[, since the lowest order quantum corrections are modified by the WZW term. 
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